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Abstract 

The plane elasticity problem for a cylindrical ly anisotropic solid 
is formulated. The form of the solution for an infinite wedge-shaped 
domain with various homogeneous boundary conditions is derived and the 
nature of the stress singularity at the vertex of the wedge is studied. 
The characteristic equations giving the stress singularity and the angu- 
lar distribution of the stresses around the vertex of the wedge are 
obtained for three standard homogeneous boundary conditions. The numer- 
ical examples show that the singular behavior of the stresses around 
the vertex of an anisotropic wedge may be significantly different from 
that of the isotropic material. Some of the results which may be of 
practical importance are that for a half plane the stress state at r=0 
may be singular and for a crack the power of stress singularity may be 
greater or less than 1/2. 


1. Introduction 

The singular behavior of the stress state in an isotropic wedge- 
shaped domain under in-plane loading conditions was studied in [1] for 
various homogeneous boundary conditions. The similar problem for the 
bending of a wedge-shaped isotropic plate was considered in [2]. The 
problem does not appear to have been studied for orthotropic materials. 
Physically the wedge problem for an orthotropic material with rectangu- 
lar axes of orthotropy is, perhaps, not very practical. However, in a 
group of materials known as the cylindrically anisotropic solids, the 
problem not only is practical but could also be important. These are 
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materials such as wood which have an axisynmetric structure. Further- 
more, since such- structured materials could be prone to brittle or 
quasi-brittle fracture, the problem of stress singularities needs to 
be studied. In this paper we consider an infinite cyl indrically aniso- 
tropic wedge-shaped domain having various combinations of homogeneous 
boundary conditions. It is assumed that in the solid under in-plane 
loading, the conditions of generalized plane stress or plane strain are 
satisfied and the vertex of the wedge coincides with the center of the 
polar axes of orthotropy. 


2. Formulation in Polar Coordinates 

Consider a cyl indrical ly anisotropic wedge-shaped domain with the 
axes of orthotropy r, 0 , and z (see the insert in Figure 1). In terms 
of the engineering elastic constants, the stress-strain relations may 
be written as 

V ° ^ “e-Vz ’ 

^9 = E7r'''’9r"r * "e ‘ ''ez °z)’ 

=z ' ('Vr ■ '■ze “e "z* ' 

^re “ ^re^^®r9 * ^6z"^e^^'^9z ’ ^zr ~ "^zr^^^zr ’ (2a-c) 

From symmetry it follows that 






(3a-c) 


In the case of generalized plane stress a_ = = 0, and follow- 

z 0 Z rz 

ing [3], if we define the effective modulus E, the effective Poisson's 
ratio V, the modulus ratio 6** and the shear parameter < by 


^ > '' ' ’'"re '’er ’ ' ^r''^9 • 

_ E , 1 ^'re ''er^ 

^ ' 2 " E^ ■ E J ’ 


(4a-d) 
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the stress-strain relations become 






re 


(5a-c) 


In the plane strain case = 0, + '>26 °e 

it may easily be shown that equations (5) are still valid provided the 
definitions given by (4) are modified as follows: 


E = IjTl 






li r. 1 * _ ^r^^“''ez'^z9^ 
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(6a-d) 


For the plane problems, expressed in polar coordinates, the equa- 
tions of equilibrium and the condition of compatibility may be written 
as 


3cr^ 
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ar 
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(7a,b) 


(8) 


Substituting from (5) into (8) the compatibility condition becomes 
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If we now introduce the stress function iJ)Cr,0) defined by 


'^r “ r3r ^ r^38“^ ’ ” 3r^ ’ ^re " 3r30 V' ’ 


0 Oa-c) 


it may be seen that the equilibrium equations (7) are satisfied and sub- 
stituting into (9) we obtain 



2k 3**<}> ^ 


1 3^ X 25^ dH 2 k 

6^r^ 30^ r 3H ' 3r30^ 


1 

" 6^r^ 



1 3^ ^ 

6^r^ 3r 


0 . 


(11) 


Equation (11) is the governing equation for a cylindrically anisotropic 
solid in polar coordinates under plane loading conditions corresponding 
to the biharmonic equation for isotropic materials. 

As in [1] and [2] we will now assume the solution to be of the 
following form: 

<l>(r,0) = r^F(0,x) . (12) 


If the boundary conditions are selected to be homogeneous, the discrete 
set obtained for X and F may be considered as the eigenvalues and the 
eigenfunctions of the problem. Substituting from (12) into (11) we find 


d*+F 

de^ 


+ a 


d2F 

diF" 


+ b F = 0 


(13) 


where 


a = 2 + 2 K ^ 


b = 5 '*X(x-1)2(x-2) - X(X-2). 


(14a, b) 


The general solution of (13) may be expressed as 

a, 0 a,0 - 0,6 -a,0 

F = Ae ^ + Be ^ + Ce + De ^ 


( 15 ) 


where 


- 4 - 



«! = [ i(-a+v'a^-4b)] ^ , «2 “ C iC-a-*^a“^-4b)] (16a, b) 

The coefficients and «2 are, in general, complex and the boundary 
conditions at 9=0 and 0=0 account for the unknown constants A, B, C, 
and D, where B is the total wedge angle (Figure 1). 

In order to apply the boundary conditions the stresses and the 
displacements must be expressed in terms of the function F. From (10), 
(12) and (15), the stress components are found to be 

, , X-2/ -ai0 -026, 

a- = X(x-l)r^ ^(Ae + Be ^ + Ce + De ^ ) , 

Tr0 = -(x-l)r^‘^(a^Ae“^® + OiBe”^® - aiCe"“^® - ozDe'"^®) , 

= r^"2[A(x+a^)e“^® + B(x+a^)e°2® + C(x+a^)e'“i® 

+ D(x+a 2 )e , (17a-c) 

The expressions for the displacements may be obtained indirectly by 
referring to' the strain-displacement and the stress-strain relations. 
First we note that from (5), (10) and (12) the strains may be obtained 
as 

r “ ^ ^ di^"^ " E ^ x(x-l)F , 

Eg = - f + ^) + f >.(X-1)F , 

= . (18a-c) 

Next, we observe that in polar coordinates in terms of the displacements 
u^ and Ug, the strains are given by 

, -!!£ , _i,n \ + 

r 3r ’ ^0 r r 30 ’ re 2 ' ar r r 39 ^ ‘ 

Now assuming that the expressions for the displacements are of the 


- 5 - 


following form 


= r'"G(e,m) , = r"’H(e,m) , 

from (19) it follows that 

- m r G , - r (G + gj) 


1 ^m-1 / L, , dG\ 

= J r (mH - H + gj) . 


(20a,b) 


(21a-c) 


Comparing (18) and (21) we find the relations to determine the coeffi- 
cient m and the functions G and H as follows: 

ra r"'-’G = Cfk 0- > - f 

r“‘'(G + = /'^C- f (XF + ^) + f x(x-l)F], 

+ ^] = - (x-1) ^ . (22a-c) 

On dimensional grounds from (22) it is clear that 

m = X - 1 . U3) 


Also, imposing the physical condition that in the finite region 0 < 0 < 3, 
0<r<rQ<« the strain energy be bounded, it may easily be shown that the 
real part of m must be positive, or 

Re(x) > 1 . (24) 


Eliminating the terms involving r (22) now gives three equations 
to determine the two unknown functions G and H (in terms of F and x). 
First, by using (15) from (22a) we find 


' E6^(x-1 ) 


2 OjO 


Ba^e 


020 2 „ 2 -^ 2 ^ . 
+ Ca^e + Do2e ) 
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a,0 

Be ^ 


+ fp 


-ai0 


— rto ft 


(25) 
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Next, from (22c), (15) and (25) it follows that 

. . 1 / 3 3 “2® 3 

H(0,X) = - E^2-(r-T)'(A^' 

1 1 “1® 
■ eTxW 2(k+v)(x- 1)] (Aaie 


-cti0 


3 -“20. 
^2^ ) 


“20 “a 10 ““2®. /__! 

+ 8026 - Ca^e - 8026 ). (26) 

Substituting now from (25) and (26) it may easily be shown that (22b) 
is identically satisfied. 


3. The Characteristic Equations for X 
(i) The case of stress-free boundaries. 


Let the wedge shown in Figure 1 have the following homogeneous 
boundary conditions: 


ag(r,0) = Tj,g(r,0) = ag(r,e) = x^g(r,3) = 0 . 

Substituting from (17) into (27) we find 
A + B + C + D = 0 , 


a^A + 028 - ajC - 028 = 0 , 


Ae 




+ Be 


“26 


+ Ce 


-otiS 


+ 8e 


-OoB 


2P _ 


= 0 


. “10 n “20 /. "“l® ““20 

ojAe ^ + 02Be ^ - aiCe - 0280 


= 0. 


(27a-d) 


(28a-d) 


The determinant of the coefficients A^(x) of the homogeneous algebraic 
equations (28) is the characteristic function of the problem. The 
characteristic equation for x may then be obtained as 


2 (“1+02)6 -(“1+02)6 

Ai(x) = -8 “i “2 - (“1-02) [e + ® ] 

2 (ai-“ 2)0 -(ai-“ 2 ) 0 , 

+ (01+02) [e +6 3 “ 0 * 


(29) 
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Let X-, (j=l,2,...) be the roots of (29) and B^/A. , C./A., and D./A. , 

J w J J J J J 

(j=l,2,...) the corresponding solution of (28) with A. being an arbitrary 

J 

nonzero constant. From (12) and (15) the eigenfunctions of the prob- 
lem may then be found as follows: 


. o, .0 B. 0210 

4 ).-(r, 0 ) = r^3(e + a ® 


S ““ij® 


4.°i 


(30) 


where and are obtained from (16) and (14) by substituting X = \y 
Technically, the general solution of the problem may then be expressed 
as 

00 

cf.(r,0) = z A.i).(r,0). (31) 

1 J J 


(ii) The case of "clamped" boundaries. 

Let now the displacements be prescribed along the boundaries of 
the wedge as follows; 

u^(r,0) = Ug(r,0) = u^(r,B) = Ug(r,B) = 0 . (32a-d) 


Substituting from (20), (23), (25), and (26) into (32) we obtain a * 
system of homogeneous algebraic equations for A, B, C, and D similar 
to (28). In this case, the characteristic equation is found to be 


(ai+a2)0 "(ai+oi2)3 

A 2 .(X) = -8pqst - (pt-qs)2[e + e ] 

(01-02)0 -(01-02)0 

+ (pt+qs)^[e + e ] = 0 


(33) 


where 

p = oi + Oi[X-vX6^(X-l ) + 26 ^(ic+v) (X-1 )^] , 

q = o| + 02 [X-vX 6 ^ (X-1 ) + 2S^(k+v) (X-1 )^] , 

S = o^ + X - vX62(x-1) , 

t = o| + X - vX62(x- 1) . (34a-d) 
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(ii1) Mixed boundary conditions. 
In this case, we assume that 


trgCr.O) = T^g(r,0) = 0 , (35a, b) 

Uj.(r,e) = Ug(r,s) = 0 . C36a,b) 

Following the procedure outlined above, the characteristic equation 
is obtained as 

(ai+a2)g 

AgCx) = -4aiqt - 4a2PS + Cai-C2) (qs-pt) [e 

-(01+02)6 (01-02)6 

+ e ] + (01+02) (qs+pt) [e 

-(01-02)6 

+ e ] = 0 . (37) 


4. Angular Distribution of Stresses Near the Vertex 

Consider, for example, the case of stress-free boundaries. 
Assuming A f* 0 and solving (28) we find 


T - T [‘(ai-aa)® ^ + (ai+az)^ - 2oi6 ^ ] , 


di 6 


a, 3 


-O06. 


A Z 


j ~ 2 [202® - (01+02)6 + (oi-02)e ^ ] , 




-a, 6 . . do3 , . “do 6 

Z = -2o2e - (01-02)6 + (01+02)6 


(38a-d) 


Substituting from (38) into (17) the stress components may be expressed 
in the following form: 

Oq = Ar^"^f(e) , = Ar^"^g(e) , = Ar^"^h(e) . (39a-c) 
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For each characteristic value X. (38) and (17) give the corresponding 

«J 

functions f., g., and h.. 

J J J 

5. Solution and Examples 

As previously stated, the roots of the characteristic equations 
(29), (33) or (37) which are physically acceptable are those having 
a real part greater than 1 (see (24)). Thus, equations (17) or (39) 
indicate that if there is an acceptable root x with 1 < Re(x) < 2, 
for a given wedge angle g and material constants, then the stresses 
at the vertex of the wedge would be singular with the power a=2-X 
(that is, a*. •> r"“). Clearly, the characteristic root of primary 

* V 

interest is the acceptable root having the minimum real part. 

From (14) and (16) it may be observed that the coefficients ai 
and depend on 5 and < only, that is, they are independent of v 
as well as E. Consequently, as seen from (29), if only stresses are 
prescribed on the boundaries of the wedge, then X would be independent 
of V and E and the singular behavior of the stresses at the vertex 
of the wedge would depend on two dimensionless material constants 
(5 and <) only. In fact* as seen from (17), the same statanent is 
valid for the stress state throughout the entire wedge. However, if 
the boundary conditions involve the displacements as in the cases of 
clamped or fixed-free boundaries considered in Section 3, then x as 
well as the stress state in the wedge would depend on v in addition 
to 5 and K. 

The properties of the materials used in the examples are given 
in Table 1. Material 1 is a mildly anisotropic material and corre- 
sponds to drawn titanium. Material 2 is strongly anisotropic and 
corresponds to graphite. Material 3 is a typical isotropic structural 
metal. Material 4 is yellow birch [4]. The table also gives the 
dimensionless constants 5, tc, and v with subscripts a and e corre- 
sponding to plane stress and plane strain conditions, respectively. 

The calculated results giving the acceptable characteristic value 
X with the minimum real part are given in Figures 1-3. From these ’ 
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Table 1. Material constants used in the examples 


Material 

1 

2 

3 

1 ■ 

4 

5 

E^Cpsi ) 

15.07x106 

40x106 


161,850 

161,850 

Eg(psi) 

20.8x106 

1.5x106 


103,750 

103,750 

'^re 

0.196 

0.187 


0.665 

0.665 

'’er 

0.271 

0.007 


0.426 

0.426 

G„(ps1) 

6.78x106 

4x106 


35,068 

35,068 

5 

a 

0.9226 

2.2724 

1 

1.1176 

1.1176 

K 

a 

1.075 

0,9321 1 

1 

1.3156 

1.3156 

V 

a 

0.2307 

0.0361 

0.3 

0.5321 

0.5321 





2.075x106 

100,000 

''zr 




0.426 

0.426 

CD 

N 




0.451 

0.451 

5 

£ 



j 

1.1187 

1.1490 

K 

£ 



1 

1.3196 

1.4288 

V 

£ 




0.5508 j 

1 

1.0463 
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results it may be observed that: 


- Depending on the material constants, for a given wedge angle B 
the power of the stress singularity 2-x at the vertex of a cylindrically 
anisotropic wedge may be stronger or weaker than that corresponding to 
an isotropic wedge. In this respect the dominant material constant 
appears to be the modulus ratio 6 = (E^Eg)^, 6 > 1 corresponding to 
singularities stronger, and 5 < 1 to weaker singularities than the 
isotropic solid. 

- In the cases of wedges with stress-free or clamped boundaries 
(i.e., the cases (i) and (ii) discussed in Section 3) the stress state 
at the vertex of an isotropic wedge is bounded for any angle B ^ it and 
is singular for B > ir. However, as seen from Figures 1 and 2, this is 
not the case for the cylindrically anisotropic materials. 

- A more important and perhaps somewhat surprising result con- 
cerns the case of the radial crack Ci*e., B = 2tt). In the isotropic 
and rectangularly orthotropic solids, the power of stress singularity 
a at the crack tip (in - r““) is always 1/2. However, as seen 
from the figures for 6 = 2 tt, a may have any value between Q and 1. 

- The leading power of singularity a = 2-X is always real for 
stress-free and clamped-clamped boundary conditions. 

- The trends of X vs. B curves are generally similar for free- 
free and clamped-clamped cases (Figures 1 and 2). However, as it may 
be seen from Figure 2, the curves giving the roots of the characteris- 
tic equations may have more than one branch in the critical range 

1 < X < 2 (Figures 2). 

The results shown in Figures 1 and 2 are obtained for the general- 
ized plane stress condition. The difference in X obtained from plane 
stress and plane strain cases in birch (one of the woods for which 
all material - constants were readily available) is found to be neglibi- 
ble. Therefore, to establish the trend, the Material ‘5 is used as an 
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example. Material 5 has the same constants as Material 4 except for 
E^. As seen from Figure 3, unlike the isotropic materials, in cylin- 
drical ly anisotropic solids there may be a difference between the singu- 
larities for wedges under plane strain and generalized plane stress 
conditions. However, in spite of the assumed exaggerated value for 
(see the properties of Materials 4 and 5 in the table), Figure 3 shows 
that the difference may not be very significant. 

Figure 4 shows an example for the ''symmetric mode" of the angular 
distribution of stresses (in the neighborhood of r=0) in a half plane 
(3=m) with stress-free boundaries under plane stress conditions. The 
curves give the results calculated from (39) for Material 2 which is 
a strongly anisotropic solid. For comparison, the isotropic results 
calculated from the same program by letting 5=1, <=1 are shown in 
Figure 5. For the symmetric mode, the exact isotropic results are 
(see equations 39, and Figure 1 with 6 = tt) 

f(0) = sin^e , g(e) = - ^ sin 26 , h(0) = cos^0 . (40) 

The comparison of Figures 4 and 5 indicates that the differences 
between the two sets of results, particularly in h(e), may be quite 
significant. 
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0 

Figure 1. Mimmum acceptable characteristic value x in a cylindrically anisotropic 
v/edge of angle 3 with stress-free boundaries under generalized plane 
stress conditions. The curves 1 to 4 correspond to the materials 1 to 4 
given in Table 1. 





Figure 2. Same as Figure 1 with fixed displacements along the 
boundaries (the clamped-clamped case).- 
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Comparison of the minimum acceptable characteristic values 
X for Material 5 with stress-free boundaries under plane 
strain and generalized plane stress conditions. 




Figure 4. Angular distribution of the stresses for Material 2 
with 3 =tt (half plane) under generalized plane stress 
condition (f(e) ~ o^, g(e) ~ h(e) ~ a^). 
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